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Abstract 

We investigate the possibility that the dimension 2 condensate 
A 2 ^ has a non zero non-perturbative value in Yang-Mills theory. We 
introduce a multiplicatively renormalisable effective potential for this 
condensate and show through two loop calculations that a non zero 
condensate is energetically favoured. 
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Recently [1] a lot of interest has arisen concerning the possibility of a 
condensate in Yang-Mills-theory of mass dimension 2. A possible candidate 
is given by the gauge invariant operator 

( min jd A x{AV) 2 ) 

1 W} 

A = 

2 V.T. 

where = in the absolute Landau gauge [ 11] and A 2 ^ = A a A a . There 
are several phenomenological reasons [2] to believe that the groundstate of 
Q.C.D. favours a non-perturbative value for this condensate different from 
zero. Theoretically, it was shown in [1] that monopole condensation in com- 
pact Q.E.D. is related to a phase transition for this condensate. In this paper 
we would like to give further theoretical evidence that the non-perturbative 
groundstate of Q.C.D. favourizes energetically a non-zero value for this con- 
densate. For this, several problems have to be solved. 

First of all, there is the question of what we mean by the non-perturbative 
value of (A 2 ). Perturbatively, this condensate is quadratically divergent and 
Borel-non-summable because of the presence of ultraviolet renormalons. If by 
non-perturbative, we mean that part of the condensate that is proportional 
to the A 2 -parameter, this is ambiguous because it depends on an arbitrary 
summation prescription for the perturbative part [3,4]. A second problem 
is how to define a renormalisable effective potential for the local composite 
operator A 2 ^. Because the composite operator is local, new divergences are 
introduced which necessitate new counterterms that spoil an energy interpre- 
tation [5] . In this letter, we will show how a unique non-perturbative value of 
the condensate (A 2 ^) can be defined. For this condensate, we will construct a 
multiplicatively renormalisable effective potential which is unique and whose 
absolute minimum gives the non-perturbative groundstate. We will calculate 
this effective potential up to two loops and show that up to this order, the 
groundstate favours a non zero value for the non-perturbative condensate 
(A 2 ). We conclude with some numerical results for the gluon condensate 
^-(F 2 U ) and some comments. 

To define the effective potential for the non-perturbative condensate (A 2 ) 
we introduce a massterm \J(A 2 ^} in the SU(N) Yang-Mills Lagrangian in 
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(1) 



the Landau gauge. This term generates new divergences in the generat- 
ing functional for connected Greens- functions W(J). There is a quadratic 
divergence linear in J corresponding to the quadratic divergence of (A 2 ). 
As we will show, this divergence drops out of the effective potential so we 
don't have to renormalise it. There is a logarithmic divergence linear in 
J corresponding with multiplicative massrenormalization which can be can- 
celled by a counterterm ^SZ 2 JA 2 l . Finally there is a logarithmic divergence 
quadratic in J which corresponds to a new divergence in the Greens function 
(A 2 {x)A 2 (y)) c when x — > y and which can be cancelled by a counterterm 
S(J 2 /2. These counterterms are sufficient to ensure a finite renormalised 
W(J). The reader might question this on the basis of the common wisdom 
that massive Yang-Mills-theory is non-renormalizable [6]. However, the mass 
term \ JA^ is added to the Lagrangian after gauge fixing. Therefore, our mas- 
sive Lagrangian is not the one of massive Yang-Mills theory. In particular, 
the vanDam-Veltman-Zakharov [7] discontinuity theorem is not valid and we 
have a smooth J — > limit. A simple power counting argument can then 
be used to show that our new counterterms renormalise the theory. We will 
discuss the problem of unitarity at the end of this letter. 

Let us now try to define a non-perturbative value of (AV). Therefore we 
consider the massive gluonpropagator G(k 2 , J) as a function of J. Suppose 
furthermore that G is a multivalued function of J. This means that if one 
starts from the perturbative groundstate at J = caracterised by a certain 
value of (A 2 ^), makes a contour in the complex J-plane around one or more 
singularities and then comes back to J = on a different Riemann sheet, 
one can end up in a non-perturbative groundstate caracterised by a different 
value of the condensate. This situation is analogous to A0 4 theory with neg- 
ative mass term where (<t>)(J) is multivalued. The role of the negative mass 
is played by the tachyon pole, generated by infrared renormalons in the A 2 ^ 
channel. What is different is that in our case, there is no spontaneous sym- 
metry breaking and that the perturbation series around the different vacua 
are identical. Then how can we make a distinction between the perturbative 
and a non-perturbative groundstate ? For that, we need a quantity which 
is zero to all orders in perturbation theory. As a candidate, we can take 
G _1 (0, 0) which because of gauge invariance, is zero to all orders in pertur- 
bation theory. Hence we can define the perturbative gluon propagator as 
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the propagator for which lim G 1 (0, J) = 0. On the perturbative Riemann 

sheet we have G(k 2 , J) = G p (k 2 , J) and the perturbative condensate is then 
defined through : 

5<4> = l/|y^ ( * 2 - )^/(^i G( * 2 - J )- G ^- )l < 2 » 

= A p + A np (J) 

The perturbative part of the condensate A p as defined is this way, is not the 
perturbative series for ^(A 2 ) summed in some arbitrary way but the value 
of \{A 2 ^) for J = on the perturbative sheet. This perturbative value is well 
defined after regularization and contains all the quadratic divergences. The 
non-perturbative condensate is only logarithmically divergent and vanishes 
with J on the perturbative sheet. 

To construct an effective action for the non-perturbative condensate 
A np ,we consider the generating functional W(J) and do a Legendre trans- 
form with respect to J. The only way that W(J) can implicitly depend on A p 
is through the linear term in J which contains the quadratic divergences. 
However, in the Legendre transform the linear terms in J cancel so the 
Legendre transform of W(J) is implicitely only a function of A np . Gauge 
invariance plays a very important role in this. Indeed, because of gauge in- 
variance, quadratic divergences cancel in self-energy subdiagrams of W(J). So 
the only possible dependence of W( J) on A p is through the overall quadratic 
divergence linear in J. As a consequence we can forget about the perturba- 
tive condensate and use a gauge invariant regularization such as dimensional 
regularization where A p is automatically zero. Introducing counterterms 
\5Z 2 JA 2 l and 5(J 2 /2 for the logarithmic divergences linear (multiplicative 
mass renormalization) and quadratic in J (vacuum energy divergences) we 
obtain a finite renormalised functional W(J) given by : 



e- w{J) = J [dA^] exp - J d D x 



T2 



\*% + \ Z2jA l ~ (C + 5 0y + C a . f + C c , 

(3) 

To ensure a homogeneous renormalization group equation we had to intro- 
duce a new independent parameter C(/-0- Defining the bare quantities 
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40 7 1/2 4 
(i — ^3 M 

Jo = -j-J 

A3 
Z 3 

M 2 = /i' e (C + 5C)J 2 (4) 

the RGE for W(J) becomes 



(5) 



where 



d 

i 2\ d Z 2 , 

V(g 2 ,0 = /^Clso,e,Co,Jo (6) 

Because of (4) and the single valued relation between \i and g 2 (n), we can 
consider ( as a function of g 2 and we have : 

d 

^C\g ,e,j ,(o = v(g 2 ,() 

= 2 l2 (g 2 )( + 5(g 2 ) (7) 

where 

S(g 2 ) = (e + 2 l2 (g 2 )-(3(g 2 )^- 2 )5( (8) 

is a finite function of g 2 . 

In defining a finite value for the energy functional W(J) we have intro- 
duced two problems. First, since we had to introduce a new parameter (, 
there is a problem of uniqueness. Secondly, for renormalisation purposes, 
we had to introduce a quadratic term in J in the Lagrangian. Naively, one 
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expects that this will ruin an energy interpretation for the effective potential 
defined via the Legendre transform. In the case of the Gross-Neveu model 
[8] , both problems were solved by one of us in [9] . Concerning the first prob- 
lem, it is possible to choose ( to be a unique meromorphic function of g 2 
such that if g 2 runs, ( will run according to (7). Indeed, the general solution 
of (7) reads 

C(9 2 ) = C P (9 2 ) + a exp ^2 (9) 
( p (g 2 ) is the particular solution of 

P(9 2 )^- 2 C(9 2 )=^(g 2 )+5(g 2 ) (10) 

which has a Laurent expansion around g 2 = : 

( P (g 2 ) = ^ + c h + c 1 h 2 g 2 + --- (11) 

where we have temporarily reintroduced the dependence on h. Note that the 
n-loop ( p will necessitate the evaluation of the (n+l)-loop renormalization 
group coefficient functions /3(g 2 ), 72(# 2 ) and S(g 2 ). If we put a = 0, we not 
only eliminate an independent parameter but the vacuum energy divergences 
become multiplicatively renormalizable : 

((9 2 ) + 5ag 2 ,e) = Z c (g 2 ,e)((g 2 ) (12) 

Since ( is now a unique function of g 2 which runs according to the RGE, the 
energy functional W(J) obeys the homogeneous RGE : 

{4, + ^-^J A 4) wij)=o (i3) 

Therefore the composite operator 

\Z2AI - Z ( (J (14) 

has a finite and multiplicatively renormalizable expectation value A# = 
and two-point function. For J = 0, A R = on the perturbative sheet while 
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A R = Z 2 A np on the non-perturbative sheet. The effective action for A R is 
defined by 

T(A R ) = W(J) - J d A xJA R (15) 

and obeys the RGE 

("I + + ^ 2 > / dixA "sh) r(A «> = (16) 

To calculate T(A R ) one can proceed in a straightforward way by calcu- 
lating W(J) and doing the inversion. This is rather cumbersome though, 
especially for space-time dependent J. A much more efficient method which 
displays explicitely the energy interpretation of T(A R ) uses a Hubbard- 
Stratonovich transformation 



[da] exp 



1 



d D x 



9 2 



(17) 



to eliminate the \Z 2 JA 2 ^ and Z^(J 2 terms from the Lagrangian. Our energy 
functional can now be written as a pathintegral over A^ and a fields 



-W(J) 



= J[dA,][da] 



exp - 



„ / . \ a J 
C{A^a) 
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d D x 



(18) 



where the a-field Lagrangian is given by 
1 



^v) 2 + £g.f. + £c.t. 



+ 



a- 



2g 2 Z c ( 



+ -fj. 



1 wo Z 2 



6/2 



1 z 2 

, oaA a a A a a + -fx e ^-(A a u A a a ) 2 (19) 
g 2 Z c ( 8 Z C C M M 



In our new expression for W(J), J appears now as a linear source term for 
the a field so that (a) = —gA R . The inversion and Legendre transform are 
therefore unnecessary and we simply have 



r(A fl ) = F 1PI (a = -gA R ) 



(20) 



which can be calculated in perturbation theory using the background field 
formalism. 
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We have obtained a new multiplicatively renormalizable Lagrangian 
£(a, Afj) which is to all orders in perturbation theory equivalent to the orig- 
inal Yang-Mills Lagrangian. If one perturbs around a — 0, one recovers 
the original perturbation series with is well known problems such a infrared 
renormalons. If one expands around a 7^ 0, one has an effective gluon mass 
which incorporates non-perturbative effects signalled by the infrared renor- 
malons. 

To see whether the groundstate favours a ^ 0, we have calculated the 
effective potential for a up to two loops. To calculate ((g 2 ) up to two loops, 
we had to calculate the R.G. functions up to three loops. The calculations 
where done in the Landau gauge in the MS scheme in D = 4 — e using the 
tensor correction method [10] which is a new method for efficient calculation 
of multiloop Feynman diagrams. We calculated W(J) up to three loops and 
found that it could be renormalised with the counterterm —5(J 2 /2 where 



+ 




)( 



665 1 6629 1 
—„ — o + 



36 



432 



(21) 



For mass renormalisation, we found 
Z 2 = 



(g 2 N c \ 3 (g 2 N c Y /53 1 
\l6n 2 J 2e l^vr 2 J \ 8 e 2 



+ 



(g 2 N c \ 

v 16vr 2 j 



144 e 3 

and anomalous dimension : 

'g 2 N c \ 35 , fg 2 N c \ 2 



5141 1 20717 1 
+ 



_ 95 r 
~ 487, 

/11713 3C(3)\ 1 



864 e 2 V 2592 



+ 



16 J e y 



;22) 



i2(g ) = 



167T 2 / 6 



16tt 2 



449 
~2A 



+ 



'g 2 N c y / 94363 9 



16tt 2 



V 864 



jgtfS)) (23) 



For the renormalisation group function of the vacuum energy, we obtained 
using (8), (21) and (23) : 



% 2 ) 



(N 2 C - 1) 
16tt 2 



-3-( 9 - 




\ 16tt 2 



(gr 

\ 16tt 2 J V 144 



N c \ / 71551 693 



16 



C(3) 



(24) 
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Finally we solved (10) with a Laurent expansion in g 2 and found up to two 
loops : 



2, W " 1) 



16tt 2 



i6tt 2 \ 9 161 (g 2 N c \ /567343 + S 
+ \16^J { 



g 2 N c l 13 52 



2539C(3) V 
35568 J 

(25) 

From (21) and (25) one can calculate Z^, so we now have all ingredients to 
calculate C(cr, A^) to two loop order. 

We can read off the effective gluon mass in lowest order from (19) and 
(25) : 

2 13 N c 

We define a' = y N^-i a so ^ nat the background field method at one loop 
gives free gluons propagating with an effective mass ra 2 ^ = g& ' . Since at 
one loop, Zq = 1 — y (f^r) \ and using the one loop value of ((g 2 ), we have 

~ l+ 13 (g*N^ 1 _ 13 161 |ViV c y 



9 fiV 2 - 1) a' 2 
U ; 13 N„ 2 



3 V16tt 2 / e 9 52 \ 16tt 2 / 



+ -Trm(-D + #V 2 ) 



(27) 



where the trace goes over color and Lorentz indices. Because there are N 2 — 1 
gluons with 3 massive polarizations in the Landau gauge, we find in het MS 
scheme : 



64tt 2 



-9 ° 



2 5 V 

t + In 

e 6 /r 



(28) 



The divergences cancel and we obtain a finite one loop effective potential : 
9 (N 2 -l)a> 2 , 3 /Ar2 ^{jgo** 



5 161 V 
-77 - -77777- + In 



6 78 



(29) 



The two loop correction has been calculated in [10] and reads : 



AV 2 (a>) = {N l-i^'f (9 2 N C 



167T 2 y 167T 2 

49359 891 C(2) 9171 



!l, n ^-lAn^ 
4 /Z 16 



A 4 



3952 32 



16 7904 



C(3) 



(30) 



9 



where s 2 = ^C7 2 (7r/3) ~ 0.2604341 • • •. 

At one loop as well as at two loops, the perturbative vacuum a' = 
is a local maximum and a lower minimum is obtained for a' ^ 0. We can 
use the RGE to sum leading logarithms and put /I 2 = go' . Introducing the 
expansion parameter y = we find a global minimum for V(a') at one loop 
for yi = 0.19251 and at two loops for y 2 = 0.14466, independent of N c . The 
corresponding coupling constants are reasonably small : for N c = 3, a ~ 0.8 
(1 loop) or a ~ 0.6 (2 loops). Through dimensional transmutation we obtain 
non- vanishing effective gluon masses. At one loop we find mi = (g<r)i ~ 
2.05A I? g ~ 485MeV for A MS = 237MeV. At two loops, we find using the 
one loop /^-function m 2 i ~ 2.59Ajg^- ~ 614MeV and using the two loop 
/3-function, m 22 ~ l-96A^r ~ 464MeV. For the non-perturbative vacuum 
energy density and for N c = 3 we find e\ ac xs — 0.335A^ at one loop while at 
two loops we find, e 2 * c m — 1, 7A^- and e 22 c — 0.567A|^. Finally we can 
calculate the gluon condensate {^-F 2 ) by making use of the trace anomaly : 

©«, = ^ (Ft? (si) 

From the anomaly we deduce for N c = 3 that the gluon condensate is related 
to the vacuum energy density as : 

(-F 2 ) = -fe vac (32) 

71 11 

Using our numerical results for e vac , in one and two loops (with one and two 
loop /5-functions) we find for the gluon condensate : 

(-F 2 )! = 0.0031Gel/ 4 

7T 

(-F 2 ) 21 = 0.0156Ge\/ 4 

7T 

(t -?2\ _ n nncon„T/4 



7T 



F 2 ) 22 = 0.0052Ge\/ 4 (33) 



Since e vac in (32) is really the energy difference between the non-perturbative 
and the perturbative groundstate, our definition of the gluon condensate is 
in fact {^F 2 ) = (^F 2 ) np — (^F 2 ) p where the suffices p and np means taking 
the J = limit on the perturbative and non-perturbative sheet respectively. 
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In this letter we have introduced a consistent definition for the non- 
perturbative value of the local composite operator A 2 and given evidence 
through two loop calculations of a multiplicatively renormalisable effective 
potential that the non-perturbative vacuum favours a non-zero value for this 
condensate. Our calculations can only be seen as qualitative indications that 
non-perturbative values for A 2 ^ can lower the energy. Other important non- 
perturbative effects such as instantons have been left out in the calculation of 
the effective potential. Since the operator A 2 ^ from which we start in equation 
(1) is gauge invariant, our results are gauge invariant. For a non zero conden- 
sate a 7^ 0, perturbative unitarity in the gluon sector is broken. It is known 
[12], that gauge invariance and perturbative unitarity should not always go 
together. However, confinement could solve this and secure non-perturbative 
unitarity in the zero color sector. 
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